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1. Introduction 

Polarizabilities represent fundamental properties of hadrons, encoding their linear response to 
externally applied fields. Experimentally, they manifest themselves, e.g., in the non-Born part of 
the low-energy Compton scattering amplitude. While the leading low-energy response is controlled 
by the static polarizabilities found in the presence of constant external fields, at subsequent orders 
of a derivative expansion, the effective hadron Hamiltonian becomes sensitive to temporal and 
spatial structures in the applied fields. The present study focuses, specifically, on weak, spatially 
constant external electric fields E applied to a neutron. When the fields are sufficiently weak, 
nonlinear effects, i.e., terms in the effective neutron Hamiltonian of higher order than quadratic in 
the external field, can be neglected, and, to zeroth order in spatial derivatives as well as second 
order in temporal derivatives of E, the neutron mass shift can be expanded as 

Am = m{E)-m{0) = --aEE^ --YeiO-{ExE) --UevE +... (1.1) 

The static electric polarizability «£ has been investigated in a number of lattice studies, cf., e.g., 
[1-6]. Here, the goal is to perform a first evaluation of the electric spin polarizability Jei within 
lattice QCD. An effort in this direction is particularly timely in view of an ongoing experimental 
program at HlyS to measure this polarizability in isolation for the first time, albeit in a proton; 
previously, only the forward and backward spin polarizabilities Jq, Jj^ of the proton [7, 8], which 
contain Yei in combination with other polarizabilities, had been accessed experimentally. 



2. External field and four-point function scheme 

Consider a spatially constant external electromagnetic gauge field A,- generating an electric 
field Ei = doAj. It enters the lattice link variables U,- as an additional phase, U,- — > exp(/fl;^A,) • Ui, 
where a denotes the lattice spacing and q the quark electric charge matrix. To evaluate only contri- 
butions to the neutron mass shift quadratic in it is sufficient to expand the relevant neutron two- 
point function in powers of E from the outset. Thus, one can perform a fully dynamical calculation 
using existing unperturbed gauge ensembles: Expanding exp(/a^A,) = 1 + iaqAi — {a^q^/2)Aj + ... 
and inserting into the lattice action decomposes the latter into an unperturbed part and an external 
field-dependent part, S = So + Sgxt, where Sext is essentially a standard j -A coupling of the external 
field to the vector current^. As a consequence, the neutron two-point function expands as 

/.r , , ^\ _ ^\m\m\\m exp(-^o) {\ -Se,a^Sl,,ll^ . .)N^{y)Na{x) 
{Np{y)N.{x)) - ;[^^][^^][^^] expMo) {l S.. + SlJ2 + . . .) ^'"'^ 

Performing the quark integrations yields the diagrammatic representation for the contributions to 
the neutron two-point function quadratic in A depicted in Fig. [l[ In the present study, only the 
connected contributions in Fig. |I] were evaluated. The four-point function scheme at which one has 
thus arrived is particularly suited to isolate the spin polarizability Jei from the other contributions 
appearing in ([Ell): As will be specified in more detail below, each vertex in Fig. [T] is associated with 
an external field insertion A = A\e\ +A2e2', however, in view of ( |1.1[ ), only contributions ~AiA2 
are of interest, which can be isolated by populating one vertex with only A\ and the other vertex 

'Additional contact terms proportional to A?, cf. [3,4], can be excluded in the present context, as seen below. 
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Figure 1: Diagrams contributing to the neutron mass shift; crosses denote j -A insertions. 

with only A2 in each diagram (and, thus, any contact terms proportional to A? can also be disre- 
garded, as noted above). This not only discards the signals from the contributions controlled by 
Ue and a^v in ( |1.1| ), but also their statistical noise^. The very specific control over the different 
contributions afforded by the present scheme thus aids in reducing the numerical uncertainties. 
Choosing the neutron state to be polarized in the 3-direction, a suitable external gauge field is 

A = {a^i+a\t, a'^ + a^t + ajt^, 0) , (2.2) 

which generates a constant Ei and an E2 varying linearly with time, as called for by ( |1.1[ ). An 
ambiguity remains in the choice of the constants a^', ^nd fixing Ei and E2 only determines 
a\ and a^- In particular, the finite spatial lattice torus is distinct from an infinite spatial domain in 
that the constant potentials a-* cannot be wholly eliminated using gauge transformations. Due to 
the torus boundary conditions, only the residual discrete invariance A,- — > A,- + 2n/q„n„L remains, 
where L is the spatial extent of the torus and qmin the smallest unit of electrical charge in the theory. 
Fractionally different correspond to different physics; namely, they represent flavor-dependent 
Bloch momenta in view of the minimal substitution pj — (p,- — qaf). Thus, despite the dependences 
on the constants aj , 0.2 and ^2 being partially related to one another via space-time symmetries, the 
neutron mass shift is influenced by additional physics distinct from the spin polarizabihty effect, 

Am = cid -{AxE) -^YeiO-{ExE) (2.3) 

even after having restricted the calculation to contributions ~ A1A2. The additional Bloch momen- 
tum dependence^ must be disentangled from the polarizability effects in order to extract Jei- To 
this end, several choices of external fields of the type ( ^^ ) will be considered and contrasted below. 



3. Extracting the neutron mass shift in the adiabatic approximation 



In the presence of external fields of the form ( |2.2[ ), the Hamiltonian of the system under con 



sideration is time-dependent. One facet of this which deserves to be noted is the connection to 



the Bloch momenta discussed above. One can, e.g., write Ai in ( |2.2[ ) in two equivalent ways: 
Ai = aj + a\t = a\{t — to), i.e., a translation of the system by the time to corresponds to the in- 
troduction of a Bloch momentum aj = —a\to. This, of course, simply expresses the fact that the 
electric field Ei accelerates or decelerates the Bloch currents. Correspondingly, the Bloch momen- 
tum dependence of the neutron mass shift will manifest itself in a characteristic time dependence 

^Note that evaluating the difference between the mass shifts obtained with neutrons polarized in two opposing 
directions does not achieve this; that eliminates the signal, but not the noise associated with the and Kev contributions. 
■^Note that a separate dependence of the form AxE can be excluded as follows: Working in Euclidean space, each 



time derivative acquires a factor / compared to Minkowski space. The terms in {2.2 ) both contain an odd number of time 
derivatives, and Am will correspondingly be extracted from the imaginary part of the neutron two-point function. By 
contrast, any putative physical AxE dependence cannot (and will explicitly be seen not to) arise in the imaginary part. 
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of the neutron two-point function, as will be seen below. 

In general, the time evolution operator for a time-dependent Hamiltonian acquires an intricate 
time dependence. Here, an adiabatic approximation will be adopted to interpret the measured 
behavior of the neutron two-point function: Since the external field can be taken to be arbitrarily 
weak, it seems plausible to consider the case that the strong dynamics equilibrate the system rapidly 
on the scale of the temporal variation of the external field. In this situation, the behavior of the two- 
point function (now projected onto zero momentum states polarized in the 3-direction) is, at times 
t sufficiently large for excited states to have decayed, 

'1 + 74, 



Gt(p = 0,0=lTr('- 



:i-/73 75)(A^(3')A^W> 



W exp 



dt'm{t') 



(3.1) 



W 



m 



(3.2) 
(3.3) 



Inserting expansions of the normalization and the mass^ in the external fields (where the superscript 
specifies the order in the external field), 

Wo + W^^^ [A{t)] + [A(0] + . . . 

mo + m(2)[A(f)] + ... 

one obtains specifically for the contribution quadratic in the external field 

Gf\p = 0,t) = Wo exp(-moO - J dt' m(^^[A{t')]^ (3.4) 

In view of this, the mass shift m(2)[A(0] can be obtained, after dividing out the unperturbed cor- 
relator G'^^\p = 0,f) = Wb exp(— /Mo?), from the temporal slope of the ratio only at a stationary 
point of the time evolution, where the time dependences of VK^^^ [A(f)] and m^^^ [A(f)] are relegated 

to quadratic order in t. Physically, this occurs when 
the quark Bloch currents cease to flow and turn 
around into the opposite direction due to the forcing 
by the external electric field. Only at such a station- 
ary point do the strong dynamics have the opportu- 
nity to form a bona fide neutron, the polarizability 
of which can thus be extracted only in the vicinity 
of that particular point in time. 



ami 


anis 




# configs 


0.01 


0.05 


351 MeV 


448 


0.05 


0.05 


759 MeV 


425 



Table 1: A^/ = 2 + 1, 20^ x 64 MILC asqtad 
ensembles with a ~ 0.124fm used in the 
present investigation. 



4. Numerical results 

The neutron two-point function was evaluated using the MILC ensembles listed in Table 
Each configuration was HYP-smeaied and chopped into two 20^* x 32 sublattices with temporal 
Dirichlet boundary conditions. Domain wall valence quarks were employed. To improve statistics, 
averages over neutron spin in the positive and negative 3-directions were taken, as well as over the 
external field ( |2.2[ ) and its rotation by n/4 around the 3-axis. Figs. ^-Q display results for the ratio 
R2{t) = Gf\p = 0,t)/G^''^\p = 0,0 for several cases of external field. 

Consider first Figs. ^ and ^, which pertain to the case of pion mass m^^ = 759 MeV. To exhibit 
the different effects at play. Fig. ^ (left) shows the result of merely using an external field of the 



^Note that the omission of a linear term m''' in the ansatz (3.3) implies that the absence of a permanent net electric 
current or electric dipole moment in an unperturbed neutron has already been put in. 
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Figure 2: Neutron two-point function ratio R2 obtained in the nij^ = 759 MeV ensemble, in the presence of 
the external fields A ~ (aj, e2{t — fo), 0) (left) and A = (aj, e2{t — tof' 0) (right), where fo = 6a. The 
neutron source is located at f = 0, electromagnetic fields are in Gaussian units. Linear and quadratic fits, 
respectively, were performed using the fit range f/a G [4, 8]. 



form A = (fli, e2{t — to), 0), which is expected to yield a time-independent effective dynamics 



isolating the o ■ {A x E) dependence in (p3|). Indeed, the correlator ratio /?2(0 exhibits linear be- 



havior, corresponding to a constant neutron mass shift determined by the slope of Rzit), cf. ([3.4|). 
Continuing with the more complex external field A = (ai, 62(1 — toy/2, 0), the expected behavior 
of the correlator ratio R2{t) is parabolic in the vicinity of the stationary point t = to, since now the 
electric field entering the d-{AxE) dependence of ( [2.3[ ) is linear in time, crossing zero a.tt = to- On 
the other hand, an additional mass shift proportional to a -{Ax E) would manifest itself as an addi- 
tional linear term in R2{t), with the effect of shifting the parabola away from t = to- Fig. ^ (right), 
depicting the corresponding numerical result, corroborates the expected behavior; there is no shift 
of the parabolic time dependence away from t = to, and therefore no additional dependence on the 
combination o ■ {A x E), as expected on general grounds, cf. footnote |3| further above. One can 
moreover verify that the quadratic coefficient of the parabola in Fig. § (right) is compatible with 
the slope in Fig. ^ (left) within numerical error, consistent with the above interpretation. 

Finally, employing the full external field A = 
{ei{t — to), e2{t — toY /I, 0), one would ex- 
pect to see cubic behavior of R2 (?) around the 
stationary point t = to, since now the time de- 
pendences of A and E in a ■ (A x E) combine 
to render the latter proportional to {t —to)^. 
If this were the only contribution to the neu- 
tron mass shift, the slope of /?2(0 at the inflec- 
tion point t = to would vanish; on the other 
hand, an additional contribution to the mass 
shift proportional to G ■ {E x E) would mani- 
fest itself as an additional linear term in R2{t), 
the slope of which could therefore be read off 
at the inflection point of the overall cubic be- 
havior. This slope, of course, determines the 
electric spin polarizability Jei of the neutron 
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Figure 3: Neutron two-point function ratio R2 ob- 
tained in the niji = 759 MeV ensemble, given the ex- 
ternal field A = (ei(f -fo), e2{t-tQf/2, 0), where 
to = 6a. The neutron source is located at f = 0, and 
electromagnetic fields are in Gaussian units. Cubic 
fit was performed using the fit range f /a e [4,8]. 
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Figure 4: Correlator ratios R2 analogous to Fig. 2, except in the nij^ = 357 MeV ensemble, with fits per- 
formed using the fit range t/a £ [3,5]. 

which is the ultimate objective of the analysis. Indeed, the corresponding numerical result for /?2(0 
shown in Fig. ^ exhibits the expected behavior, including a nonvanishing slope at the stationary 
point t = to- From this slope, one extracts the value Yei = 0.0057(8) • 10^^ fm'^, from connected 
contributions only, at the pion mass niji = 759 MeV. 

Turning to Figs. ^ and ^ which present analogous results for the lighter pion mass niji = 
357 MeV, the principal difference lies in the substantially larger statistical uncertainties, which 
preclude a clear identification of the functional form of the neutron correlator beyond the time 
f = 5a in the plots. As a result, fits in a symmetric time range around the stationary point t = 6a, 
as employed in the m^c = 759 MeV case, are not viable at the present level of statistics. Instead, 
pending improvement of the numerical accuracy, a fit in the time range t/a € [3,5] was used to 
identify the preliminary value Yei = 0.031 (21) • 10^"* fm^, from connected contributions only, at the 
pion mass nij^ = 357 MeV. Nevertheless, the data depicted in Figs. ^ and ^ display behavior entirely 
analogous to the one seen in the heavier pion mass case for small times, where the numerical 
fluctuations are under control. 

Finally, Fig. ^ relates the two data points for the spin polarizability Jei extracted above to 
partially quenched chiral perturbation theory [9], adjusted to match the present connected diagram 
calculation. The corresponding expression depends on a number of low energy constants, all but 
one of which are comparatively well determined from either experiment or previous lattice calcula- 
tions on the same ensembles. The lone poorly 
determined parameter, gi, has been adjusted 
such as to generate the two curves in Fig. ^ 
which pass through the upper and lower ends 
of the error bar of the Yei measurement at 
nijc = 357 MeV This yields a rather stringent 
estimate of ^1 = -0.15 ^o:oIto:?o' where the 
first uncertainty quantifies the spread depicted 
in Fig. ^ and the second one was obtained by 
varying the other low energy constants within 
reasonable bounds. Note that errors due to the 
truncation of the chiral expansion were not es- 
timated. 



c\j 




Figure 5: Correlator ratio R2 analogous to Fig. 3, 
except in the nij^ = 357 MeV ensemble, with fit 
performed using the fit range f /a e [3,5]. 
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5. Summary 



The present investigation provides first lattice QCD results for the electric spin polarizability 
"^i of the neutron, albeit including coimected diagrams only. The extraction of this quantity was 
faciUtated by the adopted four-point function approach, which allows for the exact eUmination of 
other contributions to the neutron mass shift in the external electric field considered, namely, ones 
associated with the static polarizability Ue and the dispersion polarizability Uev- In addition, the 
effects of constant external gauge fields, corresponding to quark Bloch momenta on a finite spatial 
torus, had to be carefully disentangled from the spin polarizability effect itself. A clear signal for 

the connected contribution to Yei was ob- 
tained at the pion mass m^j = 759 MeV; also 
a preliminary extraction at nij^ = 357 MeV 
proved possible, although a refinement of the 
statistical accuracy is desirable. The most 
striking feature of the connected contribution 
to Yei obtained here is its smallness compared 
to the result expected for full QCD from chi- 
ral perturbation theory [9], namely, a negative 
Yei larger in modulus by roughly two orders 
of magnitude. This suggests that the electric 
spin polarizability Yei of the neutron is domi- 
nated by the disconnected contributions. 
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Figure 6: Measured lattice data in relation to par- 
tially quenched chiral perturbation theory, cf . text. 
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